Stability of equilibria for the SO (4) free rigid body 
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Abstract 

It is shown that for the generalized rigid body certain Cartan subalgebras (called of coordinate 
type) of so{n) are equilibrium points for the rigid body dynamics. In the case of so(4) there are 
three coordinate type Cartan subalgebras which on a regular adjoint orbit give three Weyl group 
orbits of equilibria. These coordinate type Cartan subalgebras are the analogues of the three axes 
of equilibria for the classical rigid body in so (3). The nonlinear stability and instability of these 
equilibria is determined. In addition to these equilibria there are others that come in curves. It is 
shown that these curves of equilibria are nonlinearly stable in the sense that the only possible drift 
in is the direction of the curve itself. 

1 Introduction 

The goal of the present work is to find the analogue of the long axis-short axis stability theorem for 
the S0(4)-free rigid body. To do this, one needs to determine first what are the analogues the usual 
three axes of equilibria in the three dimensional case. It will be shown that they are replaced by special 
Cartan subalgebras that we shall call coordinate type Cartan subalgebras. For the general case of SO(n) 
it is proved that these coordinate type Cartan subalgebras are equilibria. 

If n = 4 then, on a regular adjoint orbit, all the Cartan type equilibria are organized in three 
Weyl group orbits. The nonlinear stability and instability for these equilibria is determined taking into 
account the symplectic geometry of the orbit and the complete integrability of the system. The results 
in this paper complete and extend some previous work of Feher and Marshall [S] and Spiegler [17] . 

In addition to the Cartan type equilibria, there are, on every regular orbit, curves of equilibria. It 
is shown that these are nonlinearly stable as a family, that is, if a solution of the S0(4)-free rigid body 
equation starts near an equilibrium on such a curve, at any later time it will stay close to this curve but 
in the direction of the curve itself it may drift. 

The implication of the topological structure of the energy-mementum level sets on bifurcation phe- 
nomena in the dynamics was extensively studied by Oshemkov jl2j . |13j . Bolsinov and Fomenko 



2 Equilibria for the generalized rigid body 

The equations of the rigid body on so(n) are given by 

M=[M,n], (2.1) 

where fl e so(n), M — ft J + Jil E so{n) with J — diag(Ai), a real constant diagonal matrix satisfying 
Xi + Xj > 0, for all i, j = 1, . . . , n, z 7^ j (see, for example, [16 ). Note that M = [rriij] and — [cuij] 
determine each other if and only if Xi + Aj > since mij = {Xi + Xj)uJij which physically means that 
the rigid body is not concentrated on a lower dimensional subspace of M" . 



It is well known and easy to verify that equations (2.1) are Hamiltonian relative to the minus Lie- 
Poisson bracket 

{F, G}{M) i Trace(M[VF(M), VG(M)]), (2.2) 
and the Hamiltonian function 

H(M) Trace(Am). (2.3) 

Here F,G,H E C°°(so(n)) and the gradient is taken relative to the Ad- invariant inner product 

-^Trace(Xr), X,Ye5o{n) (2.4) 

which identifies (so(n))* with so{n). This means that F — {F,H\ for all F g C°°(so(4)), where 



{•,•} is given by ( |2.2[ ) and H by (2.3), if and only if (2.1) holds. Note that the linear isomorphism 
X S so(n) I— >■ XJ + JX e so(7i) is self-adjoint relative to the inner product (2.4| and thus VH{M) = fJ. 

Let Eij be the constant antisymmetric matrix with 1 on line i and column j when i < j, that is, 
the (fc, Z)-entry of Eij equals {Eij)ki = 5ki5ij — SkjSu. Then {Eij \ i < j} is a basis for the Lie algebra 
5o{n). We have 

{EijEks)ab ~ ^aiSjk^bs — ^aj^ik^bs — ^aiSjs^bk + Saj^is^bk 

and hence Efj is the diagonal matrix whose only non-zero entries —1 occur on the ith and jth place. In 
addition, if i < j and A; < s, we get 

[Eij,Eks\ — SjkEis + SisEjk ~ SikEjs ~ SjsEik (2.5) 

where Erp :— —Epr, if r > p. 
Since 

[M, vt] = [nj + .m, n] = njn + - n^j - nm = [j, n^]. 

we see that M is an equilibrium if and only if [J, ft^] = 0. Assuming that all Xi are distinct, this 
condition is equivalent to the statement that fl^ is a diagonal matrix. 

Theorem 2.1. Let f) C So{n) be a Cartan subalgebra whose basis is a subset of {Eij \ i < j}. Then any 



element of \] is an equilibrium point of the rigid body equations (2.1) 



Proof. We have to prove that for any M e [} the matrix 51^ is diagonal. Since the linear isomorphism 
-H- M is given by a diagonal matrix in the basis {Eij \ i < j} of so{n) it follows that M G f) if and 
only if fi e t). 

fe 

So let e () with fl — J2 cesEs, where k := [n/2] = dimf) and {Ei, ...,Ek} C {Eij \ i < j} is the 

6-1 

basis of (]. Then 

/ k \ 2 k 

= II "^^^ = + H oiiap{EiEp + EpEi). (2.6) 

\s=l / s=l i^p 

Since f) is a Cartan subalgebra we have [Ei,Ep\ = which is equivalent to EiEp — EpEi for any 
l,p e {!,..., k}. Then [EiEpf = E^El = {-^ifEpEi = EiEp. Consequently, the matrix EiEp is 
symmetric. Since Ei, Ep £ {Eij \ i < j}, we distinguish the following cases for I ^ p: 
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(a) El = Eij,Ep = Ejs,i < j < s, in which case the product EiEp is not symmetric because the 
{i, s)-entry equals 1 and the (s, z)-entry vanishes. So this case cannot occur. 

(b) El = Eij,Ep = Esj,i < j,s < j,i 7^ s. Then, EijEsj is not symmetric because the (i,s)-entry 
equals —1 and the (s,i)-entry vanishes. So this case cannot occur. 

(c) El = Eij,Ep = Eis,i < j,i < s,j ^ s. Then EijEis is not symmetric because the (j, s)-entry 
equals —1 and the (s,j)-entry vanishes. So this case cannot occur. 

(d) El = Eij,Ep = Eks,i < j,k < s, {i,j} n {fc, s} = 0. In this case EiEp = Or, 



Thus, the only possible case in (2.6) is (d) which implies that 

s=l 

which is a diagonal matrix. ■ 

We shall call a Cartan subalgebra as in Theorem[l]a coordinate type Cartan subalgebra. The dynamics 
of ( |2.1[ ) leaves the adjoint orbits of SO(n) invariant. Since the intersection of a regular orbit (that is, 
one passing through a regular semisimple element of so(n)) with a Cartan subalgebra is a Weyl group 
orbit (see, e.g. [8]), we conclude that the union of the Weyl group orbits determined by the coordinate 



type Cartan subalgebras of so(n) are equilibria of (2.1) 



3 The adjoint orbits of so (4) 

The Lie algebra of the compact subgroup S0(4) ^ {A e 0((4,M) | A*- A = /4,det(^) = 1} of the special 
linear Lie group SL(4,IR) is so(4). In this section we present the geometry of the (co)adjoint orbits of 
S0(4) inso(4). 

We choose as basis of so (4) the matrices 



E, ^ 



E4 = 











" 







-1 





1 






























1 ' 


























-1 












; E2 = 



E. = 



0010 
0000 
-1000 
0000 

0000 
0001 
0000 
0-100 



Ek — 



0-100 

10 









1 

0-10 



and hence we represent so (4) as 



so(4) 



M 






-X3 


X2 


yi 







-Xi 


2/2 




Xi 





2/3 


-yi 


-y2 


-2/3 






xi,X2,X3,yi,y2,y3 G 



> . 



(3.1) 



This choice of basis was made for computational convenience as we shall see below. Note that Ei = 
-E23, E2 = Ei3, E3 = -£^12, Ei = Eu, E5 = E2i, Ee = E34,. From ([2^, it follows that the 
multiplication for this basis of so (4) is given by the following table (the convention is to calculate [row, 
column] ) : 



[•,•] 


El 


E2 


E3 


Ei 


E5 


Ee 


El 





E3 


-E2 





Ee 


-E5 


E2 


—E3 





El 


~Ee 





Ei 


E3 


E2 


-El 





Eb 


-Ei 





E4 





Eq 


-E5 





E3 


-E2 


E5 


— Eq 





Ei 


-E3 





El 


Ee 


E5 


-E4 





E2 


-El 
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In the basis {Ei, . . . , Eq}, the matrix of the Lie-Poisson structure (2.2) is 






-2:3 


X2 





-2/3 


2/2 







-Xl 


2/3 





-2/1 




Xl 





-2/2 


2/1 








-2/3 


2/2 





-a;3 


a;2 







-2/1 


2^3 







-2/2 


2/1 





-X2 








(3.2) 



Since rank so (4) — 2, there are two functionally independent Casimir functions which are given respec- 
tively by 



Ci(M) := 



1 



Trace (M^ 



and 



C2{M) 



Pi(M) 



3 

i=l 



Thus the generic adjoint orbits are the level sets 

Orb,,,, (A/) = (Ci X C2)-i(ci,C2), 



(ci,c2) e 



Note that if M ^ 0, then dCj{M) ^ for j = 1, 2. 

The Lie algebra so (4) = so (3) x so (3) is of type Ai x Ai and, consequently, the positive Weyl 
chamber, which is the moduli space of (co)adjoint orbits, is isomorphic to the positive quadrant in 
M^. In the basis of so (4) that we have chosen above, the positive Weyl chamber is given by the set 

{(C1,C2) GM^ I > l^^l}. 

To characterize the adjoint orbits of S0(4) it is convenient to split so(4) = VuffiV^, where Vu {u := 
(mi,M2,M3) I Ui ■= Xi +yi,i = 1,2,3} = M^, K := {v := (wi,V2,W3) I := Xi - y.,,i = 1,2,3} = M^. 
Instead of the independent Casimir functions Ci, C2 we consider the following two independent Casimir 
functions 

I?u(M) 2Ci{M) + 2C2(M) = Ijuf, D^{M) := 2Ci{M) - 2C2(M) = ||v||2. 

Note that 

Orbc,c,(M) = (Ci X C2)-^(ci,C2) = (Du X Dv)-^(di,d2), where di = 2ci + 2c2, = 2ci - 2c2. 

These considerations yield the following characterization of the SO(4)-adjoint orbits. 

Theorem 3.1. Denote by the sphere in M.^ of radius r. If ci > and ci > \c2\, then the adjoint 
orbit OnhcTcAAI) equals S'^^ — ^-^ — x S'^^ ^ — , where S'^p, — ^ — C Vu, S'^p, ^ — C V^, and hence it is 

' ^ V2ci+2C2 V2ci-2c2' v'2ci+2c2 \/2ci-2c2 ^' 

regular. If ci = \c2\ > 0, then the adjoint orbit OrbciC2(-^^) either 82^^ x {0}, with S^^^ C Vu, or 
{0} X S^^y^, with S^^y^ C Vv, and so it is singular. If ci = C2 — 0, then the adjoint orbit OrbciC2 *s the 
origin of so {-i) and so it is singular. 

In all that follows we shall denote by Orbci;c2 the regular adjoint orbit OrbcjC2i where Ci > and 
ci > |c2|, which is equivalent to di, 1^2 > 0. 

Using the Lie bracket table in the chosen basis given above, it is immediately seen that the coordinate 
type Cartan subalgebras of so (4) are ti,t2,t3, where 



ii ■.^spw{Ei,Ei) = <^ Mil 












b 








—a 








a 








-b 












a,b £ 
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t2 := spa-n{E2,Er^) 



-f2 



M, 



Mf. 



a 

6 

-a 
0-600 

-a 

a 

5 

-& 



a,b e 



a,b e 



The intersection of a regular adjoint orbit and a coordinate type Cartan subalgebra has four elements 
which represents a Weyl group orbit. Thus we expect at least twelve equilibria for the rigid body 
equations (2.1) in the case of so(4). Specifically, we have the following result. 

Theorem 3.2. The following equalities hold: 

(i) ti n Orb,,;,, - {M,\fc,Ml, _„Mi,,Ml, .,}, 

(ii) t^norb,,;,, = {mI„m\_„mI^,m\_^], 

(iii) tan Orb,,;,, = {A^,%,Mi,_„M3,,M3^_^}, 
where 



1 

1 



{Vci +C2+Vci - C2) 



(3.3) 



& = ^ (\/ci + C2 - - C2) . 



Proof. Let M^^ e ti n Orb,,;,,. Then Af^^^ e ti, 2ci = 2Ci(Afi_^) = + /J^, and C2 = ^2(^/1 ,3) 
a/3. This system of equations has the solutions 



where 



(a, (3) e {(a, 5), (-a, -6), (5, a), (-6 - a)}, 



a = (Vci + C2 + \/ci - C2) 



1 

71 



(Vci + C2 - - C2) . 



Similar arguments with obvious modifications prove assertions (ii) and (iii). ■ 

The intersections ti n Orb,,;,,, t2 n Orb,,;,,, is n Orb,,;,, are Weyl group orbits. 



4 The 5o(4)-rigid body 

We shall work from now on with a generic so(4)-rigid body, that is, Aj + Xj > for i ^ j and all 
\i are distinct. The equations of motion are hence M — [A/, fi], where M = Jil + ilJ, H. S so(4), 
J = diag{Ai, A2, A3, A4}, Ai, A2, A3, A4 € M. The relationship between Q = [ujij] e so(4) and the matrix 
M S so (4) in the representation ( |3.1| is hence given by 

(A3 + A2)w32 = xi (Ai + A3)a;i3 = X2 (A2 + Ai)w2i = .T3 

(Ai + A4)wi4 = yi (A2 + A4)W24 = J/2 (A3 + A4)W34 = TJs 
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and thus the equations of motion (2.1 ) are equivalent for n = 4 to the system 



ii = 


/ 1 


1 


X2X^ -f 


/ 1 


1 


1 2/22/3 




Al -I-A3 


y A3+A4 


A2 +A4 


^2 — 


( 1 


1 ^ 


XiX^ -f 


( 1 


1 


) 2/12/3 


V A2+A3 


Ai+A2^ 


V A1+A4 


A3+A4^ 


±3 = 


( 1 


1 ^ 


X1X2 4 


( 1 


1 


) 2/12/2 


V A1+A3 


A2+A3^ 


y A2+A4 


Ai+A4^ 


2)1 = 


( 1 


1 ^ 


X2y3 + 


( 1 


1 ^ 


a;32/2 


A3+A4 


A1+A3 ; 


V A1+A2 


A2 + A4 ^ 


m = 


1 


1 ^ 


a;iy3 + 


( 1 


1 ^ 


2^3 2/1 


A2 + A3 


A3+A4 ) 


V A1+A4 


A1+A2 ) 


2)3 = 


( 1 


1 ^ 


2:12/2 + 


/ 1 


1 ^ 


a;22/i 


. A2 + A4 


A2+A3 ; 


V A1+A3 


A1+A4 ) 



(4.1) 



The Hamiltonian (2.3) has in this case the expression 



H{M) = -iTrace(M17) 

- - ( ^ + ^ 

2 \A2 + A3 ^ Ai + A3 



1 



Ai + A2 



Ai + A4 



A2 + A, 



-vl 



A3 + A. 



2 

-2/3 



The Hamiltonian nature of system (4.1 ) can be checked in this case directly, writing (ii, ±2, ±3, 1/1, 2/2, 2/3)^ — 
r_(ViJ)^, where the Poisson structure r_ is given by (3.2) and 



Xl 



X2 



XZ 



2/1 



2/2 



2/3 



A2 + A3 Ai + A3 Ai + A2 Ai + A4 A2 + A4 A3 + A4 



Theorem 4.1. // £ denotes the set of the equilibrium points of (4.1 1, then £ = ti U t2 U ts U s+ U s_, 

where s± are the three dimensional vector subspaces given by 



s± := spang 



1 



1 



Ai + A4 A2 + A3 

the right 
diagonal) yields the system 



1 



A2 + A4 



-E2± 



1 



Ai + A3 



-E, 



1 



A3 + A4 



-E.± 



1 



-Ef 



Ai + A2 

Proof. Setting the right hand side of system (4.1) equal to zero (which is equivalent to Vi^ being 



f (A2- 


FA4)(A3- 


\- \i)x2Xz ~ 


{Xl- 


f A2)(Ai ^ 


'r A3)?/22/3 


^ 


(Ai- 


FA4)(A3- 


\- \i)xiXz - 


(Xl- 


f A2)(A2 - 


t- A3)?/i2/3 


= 


(Ai- 


1- A4)(A2 - 


\- \i)xiX2 - 


(Ai- 


f A3)(A2 ' 


f A3)?/i2/2 


= 


(Ai- 


1- A2)(A2 - 


\- XA)x2y3 - 


(Ai^ 


hA3)(A3H 


h Xi)x3y2 


= 


(Ai- 


FA2)(Ai- 


F A4)2:iy3 - 


(A2^ 


hA3)(A3H 


h Xi)x3yi 


= 


. (Ai- 


1- A3)(Ai - 


1- \4:)xiy2 ~ 


(A2^ 


hA3)(A2H 


h X4)x2yi 


= 



(4.2) 



Assume that xi ^ and solve for X2, 0:3, j/2, 2/3- 

(Al + A3)(A2 + A3) J/iJ/2 



X2 = 



X3 = 

2/2 = 

2/3 = 



(Al + A4)(A2 + A4) Xl 

(Al + A2)(A2 + A3) ^1^3 

(Ai + A4)(A3 + A4) Xl 

(A2 + A3)(A2 + A4) X2yi _ (A2 + A3)^ y1y2 

(Al + A3)(Ai + A4) Xl ^(Ai + A4)2 xj 

(A2 + A3) (A3 + A4) a;3yi _ (A2 + A3)^ yly3 



(Al + A2)(Ai +A4) Xl (Ai + A4)2 xj 

If j/2 = 2/3^0 then the last two equations above hold. The first two imply 2^2 = 2:3 = 0. Therefore, 
all equations in the system (4.2) hold for any yi. Obviously the system (4.2) is satisfied even if xi — 0. 
This shows that any element of ii span(i?i, E4) is a solution of the system (4.2 ). 
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Assume that at least one of 1/2, Us does not vanish. For example, if 2/2 7^ then the third equation 
in the system above yields 

A2 + A3 

-2/1 



Xi — ± 



and thus the first equation becomes 



, Ai + A3 
X2 = ±T — r^2/2- 



A2 + A4 

If also ?/3 7^ 0, then the last equation gives the same value for xi and so the second equation implies 

_ _l_ Ai + A2 
A3 + A4 

It is clear that one can let 1/2 and be zero in these relations which shows that the solution of the 
system is in this case given by the vector 

X1E1+X2E2 + X3E3 + yiE^ + 2/2 -B5 + ysEe 

A2 ~\~ A3 Ai -f- A3 Xi -\- A2 

= ±^ —yiEi ± - — — — y2£'2 ± T — -^ysEs + yiE^ + y2E^ + y^Ee 



'X1 + X4 
■ ±2/1 (A2 + A3) 

± ysiXi + X2) 



X2 + A4 
1 



A3 + A4 



Ai + A4 
1 



El ± 



A2 + A3 



1 



Ei ±y2(Ai + A3) 



-Ee 



E2 ± "T ; — ^~-^5 



A2 + A. 



Ai + A3 



^ A3 + A4 Ai + A2 

for 2/1 1 2/2 , 2/3 G ^ arbitrary. This shows that any vector in the three dimensional subspaces s± is a 
solution of the system (4.2). 

One repeats this argument for the pairs (yi, 2/3) and (j/i, 2/2) and concludes that the solutions of the 
system (4.2 1 are in the sets t2 Us_|_ Us_ and 13 Us+ Us_, respectively. Thus any solution of the system 
( |4!2l ) is in the set ti U t2 U t3 U s+ U s_. 

The converse is an easy verification: one checks that any element of ti U t2 U t3 U s+ U s_ is a solution 
of the system (4.2 1. ■ 

Note that s± are not Lie subalgebras of so(4) and that s± H ti 7^ 0, i = 1,2,3. Let us compute 
5_l_ n OrbciC2 for ci > and ci > \c2\. For an arbitrary element of s±, we have 



Xi 



0.1 



2/1 =± 



Ai +A4' 



X2 



0-2 



A2 + A3 



2/2 = ± 



A2 + A4 ' 

02 



X3 



03 



A1+A3' 



2/3 = ± 



A3 + A4 

03 



(4.3) 



Ai + A2 



Note that if — aj ^ 0, i j, i, j e {1, 2, 3}, then the equilibrium lies in tfc, where k G {1, 2, 3} \ {i, j}. 
Thus the equilibria in s± that are not in ti U t2 U t3 must have at least two of 01,02,03 different from 
zero. From (4.3) we deduce 



< ci = Ci{M) = 



\i=i i=i / 



2 V(Al + A4)2 (A2 + A3)V 2 V(A2 + A4)2 (A1 + A3) 



2 \{X3 + Xiy (Ai+A2)2 



(4.4) 



and 



C2 



= C2(M) =^a:,2/« =± 



1=1 



02 



(Ai + A4)(A2 + A3) (A2 + A4)(Ai + A3) (A3 + A4)(Ai + A2) ^ 



(4.5) 
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which shows that ci > |c2| > 0. In the generic case when at least one of the inequahties 

Ai + A4 7^ A2 + A3, A2 + A4 7^ Ai + A3, A3 + A4 7^ Ai + A2 



(4.6) 



hold, we have ci > |c2|. If, in addition, at least one of the 7^ 0, that is, the equilibrium is not 
at the origin, then |c2| > 0. Assume that we deal only with such generic rigid bodies in so(4); then 
Cl > |c2| > 0, that is, all equilibria on s± necessarily lie on a regular adjoint orbit Orbci;c2 (see Theorem 



3.1). These equilibria are not isolated. To describe them, express conditions (4.4 1 and (4.5) in the new 



variables 

61 := 



0-2 



as 



V(Ai +A4)(A2 + A3) ' ^' V(-^2 + A4)(Ai + A3) ' 

to get intersections of ellipsoids with spheres in Orbci:c2 ns± given by 
bf fX2 + X3 , Ai + A4\ , 52 fXi+Xs , A2 + A4 



V(^3 + A4)(Ai + A2) ' 



< Cl = 



A2 + A3 
Ai + A4 



Ai + A4 
A2 + A3 



A2 + A4 Ai + Aa 



Ai + A2 
A3 + A4 



A3 + A4 
Ai + A2 



C2 = ± {bj +bl+ bl) 



This shows that if C2 > then s_ n Orbci;c2 
considerations above prove the following. 



and that if C2 < 0, then s+ H Orbci;c2 = 



(4.7) 
(4.8) 
The 



Corollary 4.2. On a generic adjoint orbit Orbci;c2J Ci > \c2\ > 0, the equilibria of (4.1 1 are given by 
the twelve points in Theorem 3.2 forming three Weil group orbits in ii, i2, is, and the subsets in s± 
described by (4.7) and (4.8). 




Fi gure li The Weyl chamber. The shaded domains represent the two disjoint connected components corresponding to 
regular orbits. On a regular orbit corresponding to the upper domain we find equlibria of type ti, t2, is, B+- On a regular 
orbit corresponding to the lower domain we find equilibria of type ti, t2, is, 0-. 

The equations of the rigid body immersed in a fluid is the case of Clebsch system. It is a Hamilton- 
Poisson system on the Lie algebra e(3)and was proved by Bobenko [3], see also [14], [15], [18], that this 



system is also a Hamilton-Poisson system on the Lie algebra so (4). Nevertheless, our rigid body (4.1) 
is different from the Clebsch system. 
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5 Nonlinear stability 



In this section we study the nonhnear stabihty of the equihbrium states £ n Orbci;c2 for the dynamics 
(4.1) on a generic adjoint orbit. 



Since the system (4.1) on a generic adjoint orbit is completely integrable ([4], [6], [9], [10]), for the 
so (4) case we have a supplementary constant of motion. Using Mishchenko's method [9], [16], we obtain 
the following additional constant of the motion for the equations (4.1) commuting with H: 



I{M) = {\l + \l)xl + (A? + \l)xl + {\l + \l)xl + {\\ + \l)vl + (A^ + \l)yl + {\l + \l)yl 
Without loss of generality, we can choose an ordering for A^'s, namely 

Ai > As > A3 > A4. 

The restriction of the dynamics M = [M, il] to the regular adjoint orbit Orbci:c2 is thus a completely 
integrable Hamiltonian system 

{0Thc,;c2,^Orh,,,,^,H\orK^,,J , (5.1) 



where cjorbcj ca ^^.e orbit symplectic structure on Orbci;c2- The Hamiltonian system (5.1) has all 
equilibria given by Corollary |4.2[ These are of two types: 



Xo {M e Orbc,;c, | d (i/|orb.,..J (A/) = 0, d (/|orb.,^.J (M) = O} 
X, {M e Orb,,;,, | d (i/|orb.,..J (A/) - 0, d (/|orb.,,.J (A/) ^ 0} 

Proposition 5.1. Xq Orbc,;,, n(ii U ts U ts) and Xi = Orbci;c2 n [(s+ Us_) \ (ti U ts U ts)]. 

Proof. The proof is a direct verification. One checks that if AI G ti Ut2 Uts then d (/lorbcj c,) (^) ~ ^ 
and that if AI G (s+ U s_) \ (ti U ts U ta) then d (/|orbc :c ) (A^) ^ 0, which proves the proposition by 
Theorem O 



This is seen in the following way. By Theorems 3.2 and|4.1[ all equilibria of (5.1 ) are formed by the 



set [Orbci;c, n(ti U t2 U ta)] U [Oiha.-c^ n(s+ Us_)], where 



Orb,,;c2n(ti UtsUta) 
= {(x,y) e X 



{xi,yi) = ±{a,b),±{b,a),X2 = 2/2 = 2^3 = 2/3 = 0} 



U {(x, y) e X M3 I (3.2, y^) = ±(a, b),±{b, a),xi= y^ = X3 ^ y^ ^ 0} 
U {(x, y) e X M3 I (3.3^ y^) ^ b),±{b, a),xi= yi = a^a = 2/2 - 0} 



= {(u,v) e > 
U{(u,v) e 

u{(u,v) e 



{ui,vi) = ±(a + 6, a - b),±{a + b,b- a),U2 = V2 = = U3 = 0} 
M? I (w2, W2) = ±(a + 6, a — fe), ±(a + 6, — a), = ui = U3 = W3 = 0} 

I (w3, U3) = ±(a + a - &), ±(a + 6, & - a), wi = wi = U2 = U2 = 0} (5.2) 



(since x + y = u and x — y = v) and 



Orbci;c2 n(s+ Us_) = {intersection of the ellipsoids (4.7) with the spheres (4.8)} 



A direct computation shows that if one represents an arbitrary tangent vector to the orbit as 5AI - 
(u X m, V X n) e Tm Orbci:c2 (A/), where m, n e M'^, then 

2d (/|orb.,,.J(A/)-<5Af 

- ((A2 + \l){ui + «i) + (A? + A2)(«i - «i), (A? + \1){U2 + V2) + {\l + A2)(^.2 - V2), 

{\l + A2)(m3 + 1-3) + {\l + A^)(u3 - v^)) • (u X m) 
+ ((A2 + Al)(ui + «i) - (A? + Al)(ui - «i), (A? + a2)(u2 + 1-2) - {\l + \i){u2 - V2), 
{Xl + \l){uz + W3) - (A3 + A^)(w3 - v^)) ■ (v X n). 
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Now one can immediately see that for any equilibrium M in (5.2) we have d (/lorb^j cj) (-^) ~ ^■ 

Next, take an equilibrium in Orbci;c2 H U s_) \ (ti U t2 U ta)]. For example, assume that C2 > 
so that Orbci;c2 ns_ = 0. Therefore by (|4.3|) we have 



Ml = Ol 



= ai 



1 



1 



Ai + A4 A2 + A3 
1 1 



Ai + A4 A2 + A3 



U2 = 02 



V2 = 0,2 



1 



1 



A2 + A4 
1 



Ai + A3 
1 



M3 = as 



1^3 = 03 



1 



1 



A3 + A4 
1 



Ai + A2 
1 



A2 + A4 Ai + A3 

and at least two of ai, 02, 03 are not zero since the equilibrium is not in ti U t2 U ts. Then 
d(/|orb.,,.J(M)-<5M 

vi + A3 + A2 + A4)' (A2 - A3) (Ai - A4: 

-mi 

1712 



A3 + A4 Ai + A2 



0203 (Ai 


+ A3 + A2 


+ A4)'(A2-A3)(Ai-A4) 


(A2 + 


A4) (Ai + 


A3) (A3 + A4) (A2 + Ai) 


aifls (Ai 


+ A3 + A2 


+ A4)'(A2-A4)(Ai-A3) 


(Ai + 


A4) (A2 + 


A3) (A3 + A4) (A2 + Ai) 


0102 (Ai 


+ A3 + A2 


+ A4)' (A3 - A4) (Ai - A2) 


(Ai + 


A4) (A2 + 


A3)(A2 + A4)(Ai+A3) 


, a2a3 (A2 


-A3)(Ai 


- A4) (Ai + A3 + A2 + A4) ( 



7713 



(A2 + A4) (Ai + A3) (A3 + A4) (A2 + Ai) 
0103 (A2 - A4) (Ai - A3) (Ai + A3 - A2 - A4) (Ai + A3 + A2 + A4) 



-ni 



"2 



"3- 



(Ai + A4) (A2 + A3) (A3 + A4) (A2 + Ai) 
_ 0102 (A3 - A4) (Ai - A2) (A2 + Ai - A3 - A4) (Ai + A3 + A2 + A4) 
(Ai + A4) (A2 + A3) (A2 + A4) (Ai + A3) 

Since at least two of 01,02,03 are non-zero, it follows that d (/|orbei:cJ (^^) 7^ 0. ■ 

Since the system (5.1) is completely integrable, we have {/|orbci ca > -^lorb^j ca } ~ 0, which implies 
that at en equilibrium M G Oihc^-c^ we get 



DX 



-flOrbej;, 



(M),DX 



(M) 



0, 



where TiX 
X 



(M) and TiXu^^^^ (M) are the derivatives of the vector fields Xj^^ 



i\or^^,.,^^[^v^) ana ua^i^,.^^^ ^jw j are tne aerivati 
at the equilibrium M. Thus DX/|q_,j^^ ^ {M),'DXh\^ 



and 



(M) are infinitesimally symplectic 



relative to the symplectic form worb^j ca (^) '^^ vector space Tm Orbcj^ca- 

An equilibrium point M € 3Cq is called non- degenerate if DXj^^^-^ (M) and DX/^i^^j^ (M) 

generate a Cartan subalgebra of the Lie algebra of infinitesimal linear transformations of the symplectic 

vector space {Tm Orbci;c2, '^Orbei;o2 (M)) 

It follows that for a non-degenerate equilibrium belonging to JCq the matrices DXj^^^^ (AI) and 

DXjj^^^^^ ^ (M) can be simultaneously conjugated to one of the following four Cartan subalgebras 

-B 
A 
S 



Type 1: 



Type 2: 



Type 3: 



-A 

-B 

A 

B B 



Type 4: 



A 








0" 













-B 










A 










B 










-A 




-B 








B 




-A 
















A 


-B 










B 


A 



(5.3) 
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where A,B eR (see, e.g., [1], Theorems 1.3 and 1.4). 

EquiUbria of type 1 are called center- center with the corresponding eigenvalues for the linearized system: 
iA,~-iA,iB,-iB. 

Equilibria of type 2 are called center-saddle with the corresponding eigenvalues for the linearized system: 
A, -A, iB, -iB. 

Equilibria of type 3 are called saddle-saddle with the corresponding eigenvalues for the linearized system: 
A,-A,B,-B. 

Equilibria of type 4 are called focus-focus with the corresponding eigenvalues for the linearized system: 
A + iB,A- iB, -A + iB, -A - iB. 

The main result of the paper is the following. 

Theorem 5.2. All equilibria in ti H Orbci:c2 o*^^ ^3 ^ Orbci:c2 '^^^ of center-center type and therefore 
nonlinearly stable. All equilibria in t2 O Orbcj:c2 '^''^ "/ center-saddle type and therefore unstable. Since 
Orbci;c2 is a generic adjoint orbit, if the initial condition is close to the given equilibrium but on a 
nearby adjoint orbit, it will stay close to it for all time. 

The proof consists of analyzing each case separately. 

The equihbria in ti. We begi n w ith the study of stability and non-degeneracy for the equilibria 
M^ f, e ti n Orbci;c2 (see TheoremE2|). To determine DX^io^b^^.^^ (^^a,b) : 1^** ^ l^"- ^e compute from 
(4.1) the linearized equations at Af^^ in the variable SM £ Tj^ji ^ Orbci;c2: 



tM 

However, since 

Tm Orbei;c2(M) = kerTA/(Ci x C2) = J 6M e so (4) 



^( 1 


1 


I A2+A3 


A1+A2 




1 




A2+A3 


= 




1 


1 


I A2+A3 


A3+A4 




1 


~ ( A2+A4 


A2+A3 




A1+A4 


A3+A4 




1 


^ A2+A4 


A1+A4 


1 


1 


A1+A4 


A1+A2 


1 


1 


A1+A3 


A1+A4 



bSy2 

b6x3 
b6x2. 



(5.4) 



{xi6xt + yiSyi) = 0, ^ {xiSyi + yiSxi) = 



i=l 



where M and SM are expressed in the coordinates {xi,yi) and {Sxi,Syi), respectively (see (3.1)), we 
conclude that in the case of the equilibrium ^ we have 

T^i Orbci;c2(M) = {SM £ so(4) | aSxi + bSyi = 0, aSyi + bSxi = 0} . 



infinitesimal symplectic matrix of the linearized equations on the tangent space to the orbit is 



This means that in (5.4) the four variations Sx2, Sx^, 5y2, Sy^ are independent and thus the 4x4 

rix of 

VA2+ 

VAi + 



(ai + 



A2+A3 



+A3 A1+A4 



2+A3 A1+A2 



_ 1 
1+A4 A1+A2 









A2 + A4 A1+A4 





1 1 

A1+A4 A3+A4 





1 ]_ 

A2 + A4 A2 + A3 



A2+A3 A3 + A4 





(5.5) 



Note that the coordinates {Sx2, Sx^, 5y2, Sy^) in are not canonical Darboux coordinates for the matrix 
[w(M^ J,)] of the orbit symplectic form a;|orbei;c2 evaluated at M^ This is why the matrix above, even 
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though it is infinitesimally symplectic relative to [a;(Af^j)], does not have the usual expression of an 
element in sp(4, M). The characteristic polynomial of this matrix is 



(A2-A4)(A3- A4) , (Ai-A3)(Ai-A2) 



(A2+A3)2 V(A3 + A4)(A2 + A4) ' (Ai 

/(Ai - A2)(A2-A4) 



(Ai 



(Ai+A4)2 V(A2+A4)(Ai+A2) 
A2)(Ai-A3)(A2-A4)(A3-A4) 



(Al + A4)2(A2 + A3)2(A3 + A4)(Ai + A3) 

(Ai-A2)(Ai-A3)(A2-A4)(A3-A4) 
(Al + A2)(Ai + A3)(A3 + A4)(A2 + A4) 



A2)(Ai +A3) 
^ (Ai-A3)(A3-A4) 
(A3 + A4)(Ai +A3) 
1 

(A3 + A4)(Ai + A2) ^ 



1 



(Al +A3)(A2 + A4) 



(A2 + A3)4 (Ai+A4)4 

The discussion of the position of the four roots in the complex plane is very complicated since the signs 
of the coefficients of and vary and depend on the relative size of the real numbers a and b which 
are arbitrary. 

Thus, we proceed in a different way. We have already seen that the linear operators TiXfj^^^^ (M^ ^) 
and DX/|q^j^ i^'^ah) ^'^^ commuting, where Af^^ is an equilibrium. We shall prove that these opera- 
tors are also linearly independent. Indeed, since the equations of motion for / are 

- y2y3) 

- VsVi) 

- yiV2) 

X3y2) 

xiys) 
X2yi), 

it follows that the linear equations given by TiXj^^^^ (Af^ j^) are 





= 2 


(Ai 


-Ai; 


(a;2a;3 


±2 


= 2 


(Ai 




(a;3a;i 


±3 


= 2 


A? 


-Ai 


(xia;2 


yi 


= 2 


AI 


-A? 


{X2y3 


2)2 


= 2 


!ai 


-Ai) 


(xsyi 


2)3 


= 2 


:ai 


-Ai) 


{xiy2 



(5.6) 





= 










= 2(Ai 


-A?) 




- bSys) 




= 2 (A? 


-Ai) 


{a6x2 


- bSy2) 




= 








tM 


= 2(AI 


-Ai) 


(bSx^ 


- aSys) 


tM 


= 2(AI 


-Ai) 


{aSy2 


~ 6(5x2), 



(5.7) 



so that the infinitesimally symplectic matrix of the linearized equations on the tangent space to the 
orbit relative to the matrix [a;(Af^j)] of the orbit symplectic form is 





2a(A?-Ai) 


-26(AI-Ai) 



2a (Ai - A^ 


26 (Al - Ai 






-26(A?-Ai) 


2a (Al - Ai) 



-26(Ai-A?) 


-2a(Al-Ai) 




(5.8) 



The matrices (5.5) and (5.8) are linearly independent. This is seen in the following way. A linear 
combination of these matrices yields a 2 x 8 homogeneous linear system and we need to show that the 
rank of its matrix is 2. It turns out that three pairs among these equations are identical. Now, using 



appropriate pairs one sees that the genericity hypothesis (4.6) implies that there always exists a 2 x 2 
minor with non- vanishing determinant. 

Therefore, the span of TiXn^^^^^ i^^ab) ^^"^ '^-^i\orb i-^^ab) forms a two dimensional Abelian 

subalgebra of the infinitesimally symplectic linear maps on ^T^,/i ^ Orbcj;c2 , "^lorbcj ca ) • want to 
show that it is a Cartan subalgebra in order to conclude that that fc is a non-degenerate equilibrium. 



12 



This is the case if and only if span^ |dX^|q^j^ (-^o &)' ^^i\oib (-^^a fc)} contains an element all 
of whose eigenvalues are distinct (see, e.g. [4], §1.8.2). 

To show the existence of such an element we begin with the study of the characteristic polynomial 



+ wiz^ + wi — 0, 



of (5.8), where 



wi ^ m{Xl - Xl){\l - Xl){Xl - Xl){Xl - Xl){a^ ~ b"^)^ > since a^b 
vi ^ Sio^ + Tib^ 

S\ = 4 ^2A2A2 — A^A^ -4^ A^ -4- A^ — '^2*^4 — '^1*^2 — '^I'^s) 
T\ = 4 ~ '^i'^2 ~^ 2A^A^ — A2A^ — '^i-^s — '^S'^l "^3) ■ 



Using (3.31 we have 



VI =ci{Si+T,) + ^cl- 4{Si - Ti) 

and 

Si+Ti= i{Xl + Xl-Xl- Xlf > 

Si-Ti= 4(A2 - + - Xl){Xl - A^ + A^ - A^) > 

which shows that vi > 0. 

The discriminant of the quadratic equation in is 

Ai = vj - 4wi = 16(A? + A^ - A^ - A^)^ {S'^a"^ + T^b^ + U[b'^) , 

where 

S'l - (A? - Xlf > 0, U[ = {Xj - A2)2 > 

= 2 (^A^A^ ~ '^i'^3 ~^ 2A2A^ -|- 2A2A^ — A^A2 ~ '^2'^4) ■ 

Furthermore, 

S[a^ + ryb^ + U[b^ = 2c\{S[ + U[) + cl(T[ - S[- U[) + 2ci^ci-cl{S[ - U[) 

= 2{cl - cl){S[ + U[) + cl{S[ +T[ + U[) + 2c,^cl~cl{S[ - U[). 
Since Ai > A2 > A3 > A4 we have 

S[ + U[ = [xl - xlf + [xl - xlf > 

s[-u[ = [xj + xj-xl- xl){xl -xl + xl-xl)>o 

S[ + U[+ T{ = {Xj + Xl-Xl- Xjf > 0. 

Thus, if xl + xl ^ xl + A3 we have Ai > (recall ci > |c2|). Since vi > 0,wi > 0, the equa- 
tion t'^ + vit + wi = has two non-zero distinct negative real roots and, therefore, the equation 
z'^ + viz'^ + wi = has two distinct pairs of purely imaginary roots different from zero. Thus 
span^ |dX/^Iq^j^^^^^ (M^ j), DX/i^^j^^^^^ (M^ is a Cartan subalgebra an it is of the first type 



(5.3). It follows that the equilibrium Af^^ is non-degenerate and nonlinearly stable because it is of 
center-center type (see [1], Theorem 1.5). 

The above computations being independent of the sign and permutation of a and 6, by an analogous 
reasoning we obtain nonlinear stability for the other three equilibria ^, M\j^ _^ in the Weyl 

orbit of Mlj^. 
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If + A| = A2 + A3 the eigenvalues of DXj^^^^^ i^ab) conjugate purely imaginary of mul- 
tiplicity two. In order to determine non-degeneracy of the equilibrium we have to find a linear 
combination DX^if^^^^ i-^ab) + '^'^■^iloib (-^ah)' "where a is a non-zero real number, that has 
distinct eigenvalues. The eigenvalues of this linear combination are the roots of the equation 



where 



with 



u[ = (Ai + A2)(Ai -I- A3)(Ai -I- A4)^(A2 + A3)*(A2 + A4)(A3 + A4) > 
w[ = (Ai - A2)(Ai - A3)(A2 - A4)(A3 - XiKX.a^ + Y,a + Z^f > 0, 

Xi ^ 4(Ai + A2)(Ai + A3)(Ai + A4)'(A2 + A3)2(A2 + A4)(A3 + A4)(a2 -b^)^0 



> 



and Yi, Zi are also expressions of Ai, A2, A3, A4, a, 5. 

The discriminant of the quadratic equation u^t"^ + v'lt + w'l =0 obtained by denoting z'^ — t is 

A[ = 4(Ai + A4)*^(A2 -I- A3)^(r2a + Z2fD, 

where 

D = 2 id- cj) {{Xl - Xlf + [Xl Xlf) + c,^cl-cl{Xl -Xl + Xl- A^)(A? - A^ + - A^) 

Y2 = -2(Ai + A2)(Ai + A3)(A2 + A4)(A3 -f A4)(Ai + A4 - A2 - A3) 
Z2 = A1A4 — A2A3. 

Note that I2 = if and only if Ai + A4 = A2 + A3, because Ai + Aj > for any i ^ j- Then 
(Ai-|-A4)2 = (Aa + Aa)^ and Af + Al = A^ + A^ imply Z2 = 0. Conversely, if Z2 = 0, then A? + A| = Ai + A| 
implies that (Ai + A4)^ ~ (A2 -I- X-^)^ , that is Ai -|- A4 = ±(A2 -I- A3). Since the solution with minus is not 
possible because Aj + Xj > for any i j, we conclude that Y2 = 0. Thus, 12 = if and only if Z2 = 0. 

However, if Y2 = 0, so Z2 = which means that A1A2 — A3A4, then we also have (Ai — A4)^ = 
(A2 — As)^ and consequently Ai — A4 = ±(A2 — A3). The solution with minus is impossible because 
Ai + A2 > A3 -I- A4 since, by hypothesis, Ai > A2 > A3 > A4. Hence we must have Ai — A4 = A2 — A3 
which together with Ai + A4 = A2 + A3 implies that Ai = A2 which is also impossible since Ai > A2. 
Therefore I2 7^ and hence A'l > if we choose a ^ —^2/12- 

Furthermore, v'l has the expression 

v[ = 2(Ai + A4)2(A2 + Xsf^cj-cliXsa' + Y^a + Z3), 

where 

X3 =2(Ai+A2)(Ai+A3)(Ai+A4)'(A2 + A3)2(A2 + A4)(A3+A4)(A?-A2 + A2-A^)(A?-A2 + A2-A^) >0. 

Since u'l > and A'l > for a 7^ —^2/1^2, choosing a G M large enough we also have v'l > 0, it;^ > 0, 
and hence the matrix DX/f i-^'^a b) + oDXj i^ab) f*^^^ distinct purely imaginary 

eigenvalues (in particular, zero is not an eigenvalue). Therefore Af^^ is a non-degenerate equilibrium 
and spaujj |dXj:/|q_,j^^ ^ (-^^a&)' ^^iIoa^ ^ (-^^a fc)} ^ Cartan subalgebra of the first type in (5.3 1. 
The equilibrium is thus of center-center type and therefore it is nonlinearly stable. 

The same holds for the other three equilibria Ml^ ^, Mlj, _^ in the Weyl group orbit of ^, 

that is, all equilibria in ti n Orbci;c2 o-tg of center-center type and therefore nonlinearly stable. 
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The equilibria in £2. We proceed as in the previous case. It is easy to see that DXj^^^-^^ ^ i^ab) 
and DXh\o,^^ ,^ i^^ab) generate a 2-dimensional subspace. The eigenvalues oiTiXj^^^^^^ ^ i^'^ab) 
the roots of the equation 



where 



^ -16(A? - Xl){X\ \l){Xl - - Xl){a^ b^f < 

"2 = 4 (- Xs^Xi^ - Aa^Ag^ - Ai2A4^ + Aa"* - Ai^Aa^ + X/ + 2 Ai^Ag^) 

+ 4 (- Xs'^Xi^ + 2 A22A4' - Ai^Aa" - Xi^X^^ + Ai^ + A3'* - Aa^Ag') b\ 

The quadratic equation + V2t + W2 = has discriminant 

A2 = - 4w2 = 16(A2 - A^ + A^ - Xlf (52a* + T2a%'' + U2b'^) , 

where 

^2 - {Xl - Xlf > 0, C/2 - (A? - Xlf > 0, 
T2 = 2 (— Ai^Aa^ — Ai^A4^ + 2 Ai^Ag^ + 2 A2^A4^ — A3^A4^ — A2^A3^) . 
Moreover, the discriminant of the quadratic expression 52 + T2a^b^ + C/2^'* is 

Ti - AS2U2 = -16(A? - XDiXj - Xl){Xl - Xl){Xl ^Xl)<0 

which implies that 820.^ + T2a^b'^ + U2b'^ > and consequently A2 > 0. Therefore, the equation t'^+V2t + 
W2 = has two non-zero distinct real roots of opposite signs and thus the equation + V2z'^ + W2 = 
has two distinct real roots and two distinct purely imaginary roots. Thus fc is a non-degenerate 

equilibrium and spaujj ^DXj^^^^^^^^^^ i^ab)' ^-^iloih^-^.^^ ^^'^a fc^} ^ Cartan subalgebra of the second 
type in (5.3). Thus, M^^^ is an unstable equilibrium of center-saddle type. 

As before, the above computations being independent of the sign and permutation of a and b, 
by an analogous reasoning we obtain non-degeneracy and instability for the other three equilibria 
^, M^i^ _^ in the Weyl orbit of A/^^ and hence all equilibria in ta H Orbcj:c2 are of center- 
saddle type and therefore unstable. 

The equilibria in 13. We proceed as in case of equilibria in ti. It is easy to see that DX^^^^ (M^ ^) 
and DX/^Iq^j^^ ^ i^^ab) generate a 2-dimensional subspace. The eigenvalues of DXj^^^^^ ^ i^'^a b) 
the roots of the equation 

z'^ + v^z^ + W3 = 0, 

where 

= 16(A? - Xl){Xl - Xl){Xl - Xl){Xl - Xl){a^ - b^ > 
V3 = 4 [(A^ - XDiXj - Xj) + [Xj - Xl){Xj - XD] a'+ 

+ 4 [(A? - A^)(A? - A^) + {Xl XDiXj Xl)] b' > 0. 

The quadratic equation + v^t + W3 = has discriminant 

A3 = w| - 4u;3 = 16(A? - A^ + A2 - Xlf {S^a^ + T^a%'' + U^b'^) , 

where 

^3 = {\l - \lf > 0, C/3 = {Xl - Xlf > 0, 
T3 = 2 [{Xl - xlKxl - xl) + {Xl - xl){xl - xl)] > 0. 

which implies that 830"^ + T^a^b^ + U^b'^ > and consequently A3 > 0. 
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Since V3 > 0, W3 > 0, the equation + vat + W3 = has two non-zero distinct negative real roots 



and therefore equation + v^z"^ + W3 = has two distinct pairs of purely imaginary roots and zero is 
not a root. Thus spaujj ^DX^^^^^^^^^ (M^ j), DXj^^^^^^^^^ {M^ is a Cartan subalgebra an it is of the 

first type in (5.3 1. It follows that the equilibrium M^^ is non-degenerate and nonlinearly stable because 
it is of center-center type (see [4], Theorem 1.5). 

The same holds for the other three equilibria Af^^ Mj^ ^, M^j, _^ in the Weyl group orbit of f^, 
that is, all equilibria in n Orbci;c2 «'~e of center- center type and therefore nonlinearly stable. 

This proves Theorem |5.2[ Next, we begin the analysis of the remaining equilibria. 



The equilibria in s±. The equilibria from the families s+ and s_ are not isolated on the adjoint 



orbits. In fact, they come in curves or points described by intersecting the ellipsoids (4.7) with the 



spheres (4.8), both families having the center at the origin. 



The hnearized equations at such an equilibrium Me G Orbci;c2 ri [{s+ U s_) \ (ti U t2 U is)] are 
d 



dt 



Sxi 



-1-^x2 
dt 



dt 
d 
dt 



:Sx3 = 



5yi 



dt^y' 



dt 



Sys 



(A2 


-A3)[-(A2- 


1- A4)a3fa2 - (A3 4 


X4)a2Sx3 ± (Ai 4 


h)a3Sy2 ± (Ai - 


f A2)a2(5y3] 






(Ai + A2)(Ai 


+ A3)(A2H 


-A4)(A3 


+ A4) 




(A3 


-Ai)[-(A3- 


h A4)ai(5x3 — (Ai 4 


A4)a3fei 


±(A2 4 


Ai)ai(5y3 ± (A2 - 


F A3)a3(5yi] 






(A2 + A3)(A2 


+ Ai)(A3H 


-A4)(Ai 


+ A4) 




(Ai 


-A2)[-(Ai- 


h A4)a2(5xi - (A2 4 


X4)aiSx2 


±(A3 4 


A2)a2(5yi ± (A3 - 


t- Xi)ai6y2] 






(A3 + Ai)(A3 


+ A2)(Ai H 


-A4)(A2 


+ A4) 




(Ai 


- A4) [±(A2 - 


h A4)a3(5x2 T (A3 4 


X4)a2Sx3 


-(Ai4 


A3)a3'5y2 + (Ai - 


f A2)a2(5y3] 






(A2 4 A4)(A3 


4 A4)(A2 H 


-Ai)(A3 


+ Ai) 




(A2 


- A4) [±(A3 - 


h X4)aiSx3 =F (Ai 4 


X4)a:i5xi 


-(A24 


Ai)ai(5y3 + (A2 - 








(A3 4 A4)(Ai 


+ A4)(A3H 


-A2)(Ai 


+ A2) 




(A3 


- A4) [±(Ai - 


h X4)a2Sxi =F (A2 4 


X4)aiSx2 


-(A3 4 


A2)a2(5yi + (A3 


f Xi)aiSy2] 






(Ai 4 A4)(A2 


+ A4)(Ai H 


- A3)(A2 


+ A3) 





(see (4.1 1, (4.3)). The characteristic equation of the associated 6x6 matrix is 

t^(fc4t^ + kiaj 4 k2al + ksaj) = 0, 

where 

k4 = (A2 + A3)' (Ai + A4)' (A3 4 A4)' (A2 4 A4)' (A3 + Ai)' (Ai + A2)' > 
/ci - 4 (A3 4 A4) (A2 + A4) (A3 4 Ai) (Ai 4 A2) (A1A4 - A2A3)2 > 
A:2 - 4 (A3 4 A4) (A2 + A3) (Ai 4 A4) (Ai 4 A2) (A3A1 - A2A4)' > 
/ca = 4 (A2 4 A4) (A2 + A3) (Ai 4 A4) (A3 4 Ai) (A2A1 - A3A4)' > 

and thus there are four zero eigenvalues; recall Ai > A2 > A3 > A4. A double zero eigenvalue is 
expected since the generic orbit is four dimensional. Restricting the linearized system to the tangent 
space to the orbit (which equals kerdCi(M) n kerdC2(Af)) yields a linear system whose eigenvalues 
are the roots of the polynomial t^{k4t^ 4 kia\ 4 k2a\ + k^a^) = 0. Therefore, the linearization of 
the integrable system (5.1) on the four dimensional adjoint orbit Orbcj;c2 at an equilibrium M^. € 



Orbci;c2 n [(s+ U s_) \ (ti U t2 U Is)] has the following eigenvalues: is a double eigenvalue and there are 
two other purely imaginary conjugate eigenvalues which can also degenerate to 0. Consequently, these 



equilibria can only be degenerate cases of type 1 or type 2 in (5.3). Thus, we cannot infer any stability 
conclusion from the linearized system. 

Note that the only time that can be a quadruple eigenvalue is when 02 = 03 = fci = 0. 
As before, we use the additional constant of motion /orb^, c, Ilorh^^.^^ that commutes with 



Orbe 



H\ 



However, by Proposition 



5.1 



d/, 



Orbc 



method used for studying the stability for the equilibria in OCq 



{Me) 7^ 0, SO we can not apply the 
Orbci;c2n(tiUt2Ut3). 
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We shall use energy methods (see [T], [7], [H], [5]). If 



1 1 

Ai + A2 + A3 + A4 ' Ai + A2 + A3 + A4 

then d{H + itiqCi + noC2)(Afe) — and the Hessian Tfi{H + itiqCi + noC2)(Me) has characteristic 
polynomial 

where 

~ (Ai + A3)(A2 + A4)(Ai + A2 + A3 + A4) ^ ° 

(Ai + A4)^ + (A2 + A3)^ 

do = > 

(Ai + A4)(A2 + A3)(Ai + A2 + A3 + A4) 

^ (Ai + A2)^ + (A3 + A4)^ 

(Ai + A2)(A3 + A4)(Ai+A2 + A3 + A4) ' 

We suppose, without lose of generality, that ai ^ 0,a2 ^ 0. The Hessian D^i?orb<:j;c2 (-^^e) — D^(i/ + 
toqCi + noC2)(Me)|Tjv,^orb,j;,2, where Tm, Orhc^-c2 = kerdCi(Me) n kerdC2(Me) has eigenvalues 
0, /3i, /32, /33 computed in a conveniently chosen basis for Tm^ Orbci;c2- Using Viete's relations for the 
characteristic polynomial of D^iJorbc^ ca (-^^e) computed in the above basis we have: 





t^{t-ai){t-a2){t- 


- "3) 


= 0, 






(Ai + A3)2 + (A2 


+ A4 






(Ai 


+ A3)(A2 + A4)(Ai + 

(Ai+A4)2 + (A2 


A2 + 
+ A3 


A3 + 


A4) 


(Ai 


+ A4)(A2 + A3)(Ai + 

(Al+A2)2 + (A3 


A2 + 
+ A4 


A3 + 


A4) 


(Ai 


+ A2)(A3 + A4)(Ai + 


A2 + 


A3 + 


A4) 



A,+A2{^y+A,[-^y 

B 



where: 



= (A? - \lf{\l - \lf [(Ai + A2)2 + (A3 + A4)2] > 0; 
A2 = (A? - \lf{\l - A2)2(Ai + > 0; 
A3 = (A? - a2)2(a2 - a2)2(A2 + X,f > 0; 

B = {Xi+ A2)'(Ai + A3)(Ai + A4)(A2 + A3)(A2 + A4)(A3 + A4)'- 

•(Ai - A2)^(Ai + A2 + A3 + A4) > 0. 

This shows that /3i,/32,/33 are all non-zero and positive since D^iJorb^j C2 (-^'^e) is positive semi-definite 
as it is a restriction of the positive semi-definite bilinear form TP{H + nioCi + noC2){Me). 
Now we can conclude the following theorem. 

Theorem 5.3. Each curve of equlibria in Orbci;c2 ri [(s+ U s_) \ (ti U t2 U ts)] is nonlinearly stable. 



That is, if a solution of (5.1) starts near an equilibrium on such a curve, at any later time it will stay 
close to the curve in Orbc^ica H [(s+ U s_) \ (ii U t2 U ts)] containing this equilibrium, but in the direction 
of this curve it may drift. Since Orbci;c2 is a generic adjoint orbit, if the perturbation is close to the 
given equilibrium but on a neighboring adjoint orbit the same situation occurs. 

Remark 5.1. One can pose the legitimate question if the statement of the theorem above could be 
strengthened in the sense that the drift in the neutral direction is impossible, at least for some equilibria. 
This would then prove the nonlinear stability of such an equilibrium on these curves of equilibria. To 
achieve this, one would have to show that D^i/jorbci C2 i-^^e) definite when restricted to the leaf Lj\/^ 
(to be defined below), which would give nonlinear stability by Arnold's method (which is proved to be 
equivalent with the other energy methods, see :2:). We shall show below that the method is inconclusive 
so we do not know which, if any, of the equilibria on these curves are nonlinearly stable. 

So let's try to apply the Arnold stability method to such an equilibrium Mg. We need to study 
the definiteness of the Hessian of the constant of the motion i/orb^ c. + Q^-^Orbc c evaluated at Mg 
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restricted to the invariant level set Ljvf^ •= -^Orbc c (-^Orb<,i.c2 (-^e)) of the dynamics (5.1 ). The conditions 



d{iforb.,,., + a/orbei,.J(^e) = 0, diJorb.,,., {M,) = 0, and d/orb.^,., (Me) (since Me G Xi) imply 
a = 0. Consequently, we have to study the definiteness of the Hessian of iforbej ca restricted to 

the tangent space at Me of Lm^- We shall prove below that this definiteness does not hold. 

Let CMe(i) be the curve of equilibria for Xhq^-^^^^,^^ with CMe(O) = Me. Then Xhq^i,^^,^^ {cM^it)) — 
and by differentiation DXhq^i,^ (0) -cm^ (0) — 0- Equivalently, using the formula of the linearization of 
a Hamiltonian vector field at a critical point on a symplectic manifold, we have A^^D^i/jorb^j ^2 (-^^e) ■ 
CAfc(O) = 0, where A is the 4x4 matrix associated to the symplectic form on the adjoint orbit Orbci:c2- 
As A is nondegenerate, we obtain that D'^iJjorbci ca (-^e) ' ^m^{0) = 0, which shows that CAf^(O) is in 
the eigendirection corresponding to the eigenvalue for D^iJ|orbci,c2 (-^^e)- 

We shall prove that X/^^^^^^^^ (Me) is coUinear with CAf^(O). Suppose not; then Xi^^^^^^^^{Me) 
is not tangent to the curve cmS^) at the point Mg. Consequently, for a small s G R, along the 
integral curve of Xj^^^^^ ^ , we find an s G K such that $^°'''°i''=2 (Me) = xi, where xi ^ cmA^)- Since 

{-fforbei;e2'^Orbei;o2} = 0, tfic flows Xho^^^^,^^ aud ^^"""^i''^^ Xj^^^^^^^^^ commutc and 

hence 

Because in a neighborhood of the only equilibria for the vector field Xn^-^y^^ are of on the curve 
CAf^(t), we conclude that X2 '■— (^xi) ^ xi. But then 



( 



-H"Orl 



<J,^o-e,.2 „ <j,^0...,^.2 „$fo-=i:=. o$f°"^--^) (Me) = Me 

- -fOrh 



Since $j -"'^i;'^2 (-^j^) ^ j^^^ because Me is an equlibrium of ( |5.1| . Therefore, $s "^'"^ (-^^e) = 2^2- How- 
ever, we also have $s ''""^'"^ (-^e) = 2:1, which contradicts the uniqueness of integral curves of ^/orbe c ■ 
Therefore ^/o.be c (-^^e) coUinear with CAfe(O) ^'^"^ ^b-us in the eigendirection corresponding to the 
eigenvalue for D^_ff lorb^^.^^ (-^'^e)- 

We have Tm^Lm^ © '^lo-ch^^-^^iMe) — lA/e Orbcj;c, and, using Darboux coordinates, it is easily 
shown that we also have Xi^^^^^ ^ _L D/orbcj cji which shows that X/^^j^^ ^ (-^e) G Tm^Lm^- 

Thus the Hessian of i?orbe c restricted to the level manifold Lm^ has a eigenvalue in the direction 

^/o...,..,(Me). 

This shows that the use of the constant of the motion / does not improve the nonlinear stability 
result of equilibrium points of type %i in Theorem |5.3[ 

Remark 5.2. The eigenvalue for the Hessian is expected since the equilibrium Mg lies on the curve 



obtained by intersecting the ellipsoids (4.5) and (4.4), both having centers at the origin; the 0-eigenspace 



is tangent to this curve of equilibria, as proved above. So the only stability we can expect is the stability 
transversal to the direction of the curve of equilibria. 
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